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ABSTRACT 


Experimental investigations were undertaken to determime the collapse 
pressures of relatively thick circular cylindrical shells subjected to 
hydrostatic pressure. A total of twenty-three fiber glass reinforced 
plastic cylinders with thickness to diameter ratios ranging from 0.05 to 
0.09 was tested to determine collapse pressures and the effective modulus 
of elasticity of the material. The length to diameter ratio was kept con- 
stant. Results were plotted non-dimensionally with P/E as the ordinate 
and h/D as the abscissa. The results were compared with available 
theoretical instability formulas and were in very good agreement with 
theory for ratios of h/D up to about 0.07. Above this value of h/D, 
experimental collapse pressures were lower than those predicted by 
theory. 

It was concluded that the von Mises equation [equation (6), Ref. 22| 
is the best instability equation for cylinders with ratios of h/D greater 
than 0.05. It was also concluded that for almost all practical cases, 
the thin shell theory is sufficiently accurate for predicting instability 
failures, 

Derivations for the following are included: (1)Expressions for 
calculation of the moduli of elasticity in the circumferential and 
longitudinal directions; (2)An expression for the shear modulus and 
for the effective modulus of elasticity of an orthotropic material; 

(3)A determinant expressing the critical pressure of simply supported, 
orthotropic, circular cylinders subjected to hydrostatic pressure; and 
(4)An expression for calculating the bending stress in a simply supported 


circular cylinder subjected to hydrostatic pressure, 
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NOTATION 


EAC 


Flexural rigidity, TTT 


Mean diameter 
Inside diameter 
Outside diameter ** 
Modulus of elasticity 

Effective modulus of elasticity for anisotropic material 


Modulus of elasticity in the longitudinal direction 


Ee, Eg Modulus of elasticity in the circumferential direction 


G 


a? 


My 


Mg 
My g 


M¢~x 


a 


Shear modulus, G= BaD 


Thickness of shell 


Non-dimensional term from thick cylinder stress formula, 
mec: 

~ R2-R? 

Non-dimensional term from thick cylinder stress formula, 


Ri + R* 
Ko = ———— = K,- | 
2 Ri - Re “ 


Effective length 


Kj 


Actual length 


Number of waves or lobes along effective length of cylinder at 
time of collapse 


Strain-pressure ratio, circumferential strain gage 

Strain-pressure ratio, longitudinal strain gage. 

Axial bending moment per unit of circumference 

Circumferential bending moment per unit of length 

Twisting moment per unit of length on element cut out by generators 


Twisting moment per unit of circumference on element cut out by 
plane perpéndicular to axis 


In the notation used in this paper, the subscript "1" refers to the 


inside of the cylinder, and the subscript "2" refers to the outside of 
the cylinder, 


ay. 





Nx 


Nx¢ 


Ngx 


Cp 
Qx 


Ry 
R2 


Number of lobes or waves in a complete circumferential belt around 
cylinder at time of collapse 


Normal force per unit of circumference in axial direction 
Normal force per unit of length in circumferential direction 
Shear force per unit length on elements cut by generators 


Shear force per unit of circumference on elements cut by planes 
perpendicular to axis 


Pressure 

Radial shear in the xz plane 
Radial shear in the yz plane 
Radial coordinate 

Mean radius 

Inside radius 

Outside radius 

Axial displacement 
Circumferential displacement 
Radial displacement 


Perpendicular distance from mid-surface of shell, positive inward 


tn 
3° D 


r 





Non-dimensional term in buckling equations, « = 


pw 


Shearing strain 


Shearing strain along i-plane in j-direction 
Normal strain in axial direction 
Normal strain in circumferential direction 
Normal strain in radial direction 
Cylindrical (angular) coordinate 
tr D 


Non-dimensional term in buckling equations, N= ™-=77— 


Poisson's ratio 





Se Normal stress in axial direction 

Og Normal stress in circumferential direction 
7.. Shear stress along i-plane in j-direction 
p Cylindrical (angular) coordinate 

XX, Change of curvature of generator 

X~g Change of curvature of circumference 


Xxg Twist of element due to Myg 
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INTRODUCTION 


The absence of experimental results and elastic instability analyses 
for thick cylinders has been pointed out by Wenk1, On the basis of this 
information and consultation with Mr. John Pulos of the Structural Mechanics 
Laboratory, David Taylor Model Basin, the authors decided to undertake an 
experimental investigation of the elastic instability of relatively thick 
cylinders subjected to external hydrostatic pressure. Thus, the primary 
purpose of this thesis evolved; which is the experimental evaluation of avail- 
able theoretical instability formulas as applied to cylinders with thickness 
to diameter ratios preater than those heretofore tested. From the results 
of the proposed experiments it was anticipated that the following question 
posed by Dr. Werk (see "Discussion" of Ref, 1) could be answered: 

"At what ratio of thickness to diameter are errors in the thin-shell 
approximation unacceptable?" 

A prime requisite prior to undertaking the thesis was the making of a 
careful survey of the literature to determine the extent of theoretical and 
experimental work, The authors conducted such a survey both prior to 
beginning the experiments and on a continuing basis throughout the period 
of the thesis work. The results of this search are presented in the section 
entitled "Literature Survey." 

It was also necessary to design and construct the experimental appa- 
ratus and the models. The materials and methods of design, together with 
the procedures employed, are indicated in the body of this parer. Since the 
models were fabricated from glass-reinforced plastic, the elastic formulas 
required in analyzing the data were modified to account for the difference 
in the values of the principal moduli of elasticity. Derivationsfor the 


expressions used are presented in the appendix. 


di 
References are listed beginning on page 67. 
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The thesis begins with a section on concepts and assumptions employed 
in the theory of thin shells. This study was made in order to gain an 


insight into the differences between a thick and a thin shell theory. 





CONCEPTS AND ASSUMPTIONS IN THE THEORY OF THIN SHELLS 
ee LT UE IAIN SHELLS 


A thin shell may be defined as one for which the ratio of the thicke 
ness of the shell to the principal radii of curvature is small compared 
to unity. The problem of classifying a shell as thick or thin therefore 
reduces to the determination of a value of h/R below which solution errors 
will be within a prescribed limit. Novosti lov has stated that for or- 
dinary technical calculations it 1s admissable to have a relative error 
of 5%, and that for thin shells this corresponds to a maximum h/ Ri = Oaeen 


: (p.354) has suggested a value of h/R = 0.10 as the dividing 


Timoshenko 
point for thin curved bars. 

There are three principal assumptions employed in the theory of thin 
shells in addition to those made for most problems in elasticity. They are 


(1) The normal stresses acting on planes parallel to the middle 
surface may be neglected in comparison with the other stresses. 


(2) The normals to the undeformed middle surface are deformed with- 
out change of length into the normals to the deformed middle 
surface. 


(3) The thickness h is very small compared to the least radii of 
curvature (i.e. h/R< 1). 


Assumptions (1) and (2) essentially reduce the problem to the case of 
plane stress. Consider the generalized matrix form of Hooke's law (see 


Pane, Be 59) 
T =AI\E, + 24N 


where 
T = the stress tensor 
2 = the strain tensor 
Ij= the trace of the strain tensor 
= the unit, 3x3, matrix 
A, A = constants 
d ae S = 
(see Ny-2v) ae 2 (/+2) 





Now at the inside surface of the cylinder (i.e. at z = h/2) the 
stress component o, is zero. Then, if the shell is thin, this compo- 
nent can be assumed to be zero throughout the thickness of the shell 
since it is small compared to the stresses o . and Cae 

Assumption (2) is equivalent to assuming that the shearing strains 
Ye, and VY, , are zero, for the definition of shearing strain is the 
change in value of an originally right angle from the unstrained state. 


The assumption is that this right angle is preserved during deformation. 


Therefore, letting o, = V,, = Yo z be zero, the tensor equation 
becomes 

Ox Oxe Oxz / OE, en 2 Yeo 0 

We So der| = A(€.+€,+€,)|0 / O} + 24142%6 Eommne 

Txz Jez O 0 O / O O Es 


Expanding the equation gives 


Km = A(€,+€o + €,) + ZE, 

7 = A (Ex + Eg + €,) + La EG 
Oz O= ACE + €5+€,) 1 Zu €2 
ao < "x6 

mr, = %2,= 9 


Eliminating €, from the first two equations by use of the third gives 


the desired stress-strain relations 


ca 
x = (€. + ace) 


/-2* 
EE 
Oe hae (€, +r &) 
= 


cre ~ Z(iaw) ‘xe 


Furthermore, assumptions (1) and (3) are utilized in reducing the 


strain components to a usable form (see Wang’, pp. 335-340). When this 





is done, the equation for Ty becomes 
a 





On = se x2 [Exe Hes oz Xe a 2 ae | 
where aK 7 our 
Exo = 9x = ore 
OL a eee Ee 
oo 2a dA a. os Hones 5 Qa. SS 


The expressions for the other stress components are similar in degree 
of complexity. 
Finally, to obtain the relations between the stress resultants 


and couples, the following expressions must be integrated. 


Ne Jo Nx ox 

M. +h Z Ug Mx he 2 Ox 2 
shel fob ee fede VE (oe he 
Nex a Te x Nx AY Txe 

eA | 2 Fey Mx) 2 Tao 


Again, assumption (3) is used and 2/Re is neglected, Then, Nyo = Noy 
Myg = Moy, and the expressions may be integrated without difficulty. 


For example, the expression for Ny becomes 
Bale Oe] ie 
imme ge mere ee? “x 98 ea 


This equation and the remaining ones obtained in a similar manner are 
known as Love's first approximation. 


It is possible to derive expressions for the stress resultants and 
couples without resorting to the use of assumptions (1) - (3). Naghdi®, 
using a variational method due to Reissner’, has made such a derivation, 

In order to gain an appreciation for the simplicity of the above expression 
and of the complexities involved in a thick shell theory, the more exact 


expression for N, is given below. 








a 
Ne =\7ege lene rtool} + foSe} + {-72> elem (+r) 


ne Pa 
a a a Z [| (eed + my Jee] f 


Presumably, these more exact relations could be introduced into the 








equilibrium equations, and, provided a solution could be effected, a 
thick shell analysis for the instability problem would become available. 
An analytical approach to a thick shell instability analysis has not 


been attempted in this thesis, however. 





LITERATURE SURVEY 


The theoretical and experimental analysis of the buckling of 
cylindrical shells subjected to hydrostatic pressure has attracted the 
attention of numerous investigators dating from experiments conducted by 
Fairbairn® in 1858. All investigators employed the theory of thin shells 
in arriving at theoretical formulas and most experiments were carried out 
on models with relatively thin walls. Tables and descriptive information 
are presented in this section to indicate the extent of previous theoretical 
and experimental investigations. Particularly helpful in the making of 
the literature survey were two reports by Nash? which contain a listing 
of 2339 papers and books relevant to experimental and theoretical work 


on shells and shell-like structures from 1828 through 1956. 


Results of Previous Experimental Investigations 


Table 1 is a chronological listing of various investigators together 
with the range of geometries involved and other pertinent data. Since the 
first analytical solution to the instability problem was not available 
until 1913, all investigations prior to this time were made principally 
to determine empirical design formulas. Those investigations made in later 
years were carried out to test the validity of analytical solutions either 
in general or as applied to a particular problem. 

With the exception of Southwell, all experimental investigators prior 
to 1929 restricted the loading on the cylinder to a radial pressure only. 
The axial thrust was eliminated by the use of a relatively stiff rod which 


connected the end plugs. The use of radial pressure only resulted from the 


§ 
i 





absence of analytical formulas based on combined loading, and from the 
fact that one of the principal applications of the results was to the 
study of boiler flues, Southwell] did use a combined loading to simplify 
the apparatus. He neglected the effect of axial loading, however, in 
comparing test results with his analytical approach. 

Most investigators after 1929 employed models subjected to hydrostatic 
(combined radial and axial) loads. An analytical solution to the problem 
of the cylinder under hydrostatic pressure was available by this time, and 
the submarine presented an immediate application. 

From columns 4,5, and 6 of Table 1 the following is noted: Most of the 
models were relatively small in diameter, as was the h/D ratio. In those 
cases where the h/D ratio was greater than 0.05 the model was either so 
long that the values of the critical pressure approached that of an infinite 
tube, or they were in the range where yielding of the material occurred 
rather than an instability failure. 

The type of support used by the various investigators is indicated 
in column 7. In general, this information was procured from sketches of 
the models as prepared for testing. A clamped type of support is indicated 
only where the author concerned specifically labelled the test condition as 
such. It is again noteworthy that no attempt was made to clamp the ends 
until iis Vole when an analytical solution became available. Also, as indi- 
cated by Cook? simply supported ends can be approached experimentally 
much more closely than fixed ends, 

Column 8 indicates the maximum collapse pressures achieved by the 
various investigators; column 9 the number of models tested; and column 
10 the variety of materials which have been utilized in the construction 
of models. 

Several general observations can be made based on the survey of 


experimental results found in the literature: 


BD 





(1) 


(2) 


(3) 


Scale effect is not important provided reasonable care is taken in 
fabricating the models. After testing 100 models representing 

the strength hull of a submarine, Windenburg!* concluded that 
scale effect was not evident. Some the earlier experiments 

have not given good results when compared to theory because 

of the inadequacy of materials and unsound fabrication techniques. 
Most of the results obtained in the last two decades have given 
very good correlation. 

Although the assumed end condition is taken as either simply 
supported or as clamped, neither condition can be precisely ob- 
tained in practice. Based on tests of 258 brass models, cornel113 
concluded that the end plugs inserted into the ends of the tubes 
offered a restraint that tended to reduce the effective length. 

He found that to get close agreement of his test results with 
theory, it was necessary to use an effective free length which 
differed from the actual free length by 4% to 1%. (Because of 

the small value of such a correction, the authors assume the 
measured free length to be the same as the effective free length 
for the tests reported in this thesis.) 

Several investigators have conducted carefully controlled experi- 
ments to determine the influence of other factors, such as eccen- 
tricity, on the collapse pressure, Sturm9 was probably the first 
to investigate the problem of out-of-roundness from an experimental 
and analytical point of view. In recent years, Cleaver+4 has 
conducted very thorough, extensive, and carefully controlled 


tests on 530 models subjected to radial pressure. He lists the 





following factors as affecting the collapse pressure: 
a. Fecentricity of bore relative to the external surface. 
b. Variations from true circular shape (out-of-roundness). 


c, Ovality of the specimen (the difference in the maximum and 
minimum diameters divided by the nominal diameter). 


d. Variation of material properties, 

From the results of his tests Sturm concludes that eccentricity of bore 
is the predominant manufacturing variable affecting collapse strength, 
but for eccentricities within the limits imposed by current specifications, 
its effects are small and can be neglected for practical purposes. 
Variations in circularity exert no measurable systematic influence, and 
provided they are within current specifications (i.e. for commercial tubes), 
they too may be neglected, 

The authors point out in passing that for a well-made, simply 
supported model the above effécts tend to cancel the effect of any 
end restraint that might be present. Therefore, neglecting these effects 
may be justified when attempting to correlate experimental and theoretical 


results. 


Results of Analytical Investigations 


Table 2 presents a summary of results of principal theoretical 
investigations. The method of approach to the problem is indicated. 

The differences in the various formulas can be accounted for by 
the inclusion or omission of certain higher order terms in the differ- 


ential equations of equilibrium, by the method of approach (i.e. equil- 


fi 





ibrium, energy, kinetic ), and by the degree of simplification of the 
final result. 
All of the more exact formulas are found to have the form 
p/® = Fy (h/D)” + Fy (h/) 
where Fy and F, are functions of 
(a) The length to diameter ratio. 


(b) Poisson's ratio. 


(c) The number of lobes in a complete cicumferential belt at the 
time of collapse. 


(d) ne number of lobes along a panoaaeon of the cylinder at the 
me of collapse. Von Mises~~ has shown this quantity to be 
unity for minimum buckling pressure. 

In order to facilitate comparison of the buckling equations with 
test results, values of Fy and F, for the model geometry and material 
under consideration in this paper have been calculated for the various 
buckling formas. These values are given in Table 3. Curves for p/E 
versus h/D are also included for these cases. (See Fig. 13). 

In addition to the formulas for the buckling pressure of cylinders 
under hydrostatic load, it is of interest to examine formulas for buck- 
ling pressure under uniform radial load only. Table 4 gives some of the 
principal radial load buckling equations. A comparison of the radial 
load formulas with the hydrostatic load counterparts indicates that the 


following relationship is approximately true 


buckling pressure under hydrostatic load = ne - 1 
Se. ee ee eS Ee 
buckling pressure under radial load n“ -1+A 


2 
Thus as A approaches zero (L/D ~~), the formas for hydrostatic 
load should yield the same results as the formulas for radial load. 
For the case where A = 0, n becomes two, and all of the more exact 
hydrostatic and radial load buckling pressure equations reduce to the 


Bresse equation (Table 4) for a tube of infinite length. 
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TEST APPARATUS, MODELS, AND TEST PROCEDURE 


Test Apparatus 


The test apparatus consists of a pressure vessel, pump, tubing, 
pressure gages, and equipment required for the taking of strain gage data. 
The arrangement and details of the apparatus are shown in Figures 1 through 
5. Specifications which were drawn up for the test apparatus are given 
below. 

Size of chamber; 5" I,D., 13" depth 
Working pressure; 10,000 psi 
Penetrations; One fitting for pressure tubing, one fitting for 
four wires from strain gages on the external 
surface of the model, one fitting for 1/4" tube 
to vent model to the atmosphere and to lead out 
internal strain gage wires. 
Pump: Hand operated pump with self-contained check valve and 
reservoir--the type commonly used with hydraulic jacks. 
Pressure gages: One 5000 psi 8" Bourdon gage graduated in in- 
crements of 50 psi, and one 15,000 psi 6" Bourdon 
gage, graduated in increments of 100 psi. 
Strain Indicator; Baldwin, Type N 

The pressure vessel was constructed in accordance with a basic design 
drawn up »y Mr. Kenneth Hom of the David Taylor Model Basin. The design 
was modified slightly by the authors, to suit the purpose required. The 
following calculations show the factors of safety, for a working pressure 
of 10,000 psi. The material is low carbon steel, the yield strength of 


which was assumed to be 30,000 psi. Formulas used are given in reference 


(Ne 


nar a ee 





t = wall thickness, inches, = +<**2 for walls 

2.125 for bottom 
d = inside diameter of chamber, inches, = 5.00 
P = chamber pressure, psi 


S = maximum stress, psi, = 30,000 
R = inside radius of chamber, inches,= 2,50 


F.S. = factor of safety based on yield strength of 30,000 psi 


Based on circumferential stress in wall: 


a Ot lf = 0,000) (2.2 = 
P= —-—— = ONG ase = 17,532 psi 
R + 0.6t 2250 4 (0.6) (2.25) ° 


Zeo,. = 17,532 = 
10,000 


ie 


Based on longitudinal stress in wall; 


mee 25t 2 5 Pee ee) De 
a R = 0.4t 2.50 = (0.4) (2.25) B4,375 pas 
~ 84,375 = g, 
F.S. 10,000 LL 
Based on stress in bottom: 
ear 2 
es St ou 630,000) (2.12 = 
eenae(0.162) ~ easy eee 
~ 23:449 = 
eA 10,000 B34 


* Ref. (39), para. UG-27(c)(1) 


*® Ref, (39), Section VIII, para. UG-34(1). 
Para. UG-36(3) states: 


"Single openings in vessels not subject to rapid fluctuations in 
pressure do not require reinforcement other than that inherent in the 
construction under the following conditions: 
eeee(b) Threaded, studded, or expanded connections in which the hole 
cut in the shell or head is not greater than two-inch pipe size." 
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The pressure vessel was proof tested to 11,000 psi on 30 August, 
1960. It was raised to this pressure three times, and strain data were 
taken on the outside surface using a circumferential strain gage. This 
test was made prior to the making of the penetrations for fittings "A" 
and "B" (Fig. 4). The strain data taken plotted linearly from zero to 
11,000 psi. 

This apparatus should not be used with a working pressure in excess 
of 10,000 psi, unless the pressure vessel is again proof tested at the 
higher pressure With strain gaze measurements taken in regions of high 
stress. In any case, it is recommended that a pressure of 15,000 psi 
not be exceeded, and that for working pressures in excess of 10,000 psi 
the pressure vessel be placed in a pit or behind a suitable barrier. 

All tubing and fittings used should be rated at the maximum working 
pressure employed. 

Persons making further use of this apparatus are warned that 
valve "G", Fig. 4, is rated at 3000 psi. It is recommended that this 
valve be replaced by one rated at 10,000 psi or higher. As a safety 
precaution in performing the tests for this thesis, a 1/8" steel plate 
was attached to the table between the overators and fittings "D" through 


"H" (Fig. 4). The plate can be seen in the photograph, Fig. 5. 
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FIGURES 


ARRANGEMGNT OF FITTINGS IN PRESSURE VESSEL 
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FIGURE 5 | a 
Photograph of Assembled Apparatus 
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Model Design 
Design of the models was based on the following considerations; 


(1) Maximum working pressure of test apparatus; 10,000 psi 

(2) Dimensions of test chamber: 5" I.D. and 13" depth 

Since the purpose of the thesis is tc evaluate the instability formulas 
for relatively thick-walled cylinders, it was necessary to design the models 
to fail by instability, and to have a relatively large thickness-—diameter 
ratio. Mr. John Pulos of the David Taylor Model Basin recommended a range 
of h/D from 0.05 to 0.10. Also, Timoshenko? (p.354) has given a value of 
h/D = 0.05 as an estimated division point between thick and thin curved 
bars. Accordingly, the range of values for h/D was selected as 0.05 to 
0.10. 

A calculation revealed that a material with a relatively low modulus 
of elasticity was required in order to cause instability failures in the 
range of h/D greater than 0.05, and, at the same time, remain within the 
upper limit of test pressure of 10,000 psi. Captain Hinners suggested the 
possible use of glass-reinforced plastic, especially in view of the current 
interest in this material because of its favorable strength-weight ratio. 

The David Taylor Model Basin had published a report8 (No. 1413) on 
the hydrostatic pressure tests of a cylindrical shell of a glass fiber 
reinforced epoxy resin fabricated by the Zenith Plastics Co. of Gardena, 
California. The circumferential modulus of this material was reported to 
be 4.8 x 10° psi, and the plot of pressure vs. strain was linear with but 
a slight departure from linearity near the collapse pressure. This material 
seemed to meet the requirements well. It had disadvantages, however, in 
that it was neither homogeneous nor isotropic, thus violating two basic 


assumptions in the usua] approach to the theory of elasticity. 








It was felt, nevertheless, that these disadvantages would be overcome 

by computing an "effective modulus of elasticity" from strain data, and 
by using models fabricated in such a way as to give maximum dispersion 

of the glass in the resin. Correspondence with the Zenith Plastics Co. 
was initiated to obtain information regarding the procurement of the 
necessary models. Zenith very generously offered to supply the models at 
no cost, 

At this point peactiicat tens were drawn ur for the model geometries. 
Using DTMB formula (9) [Formula No. 5, Table 2] , a set of curves was 
drawn for various values of L/D, with buckling pressure as the ordinate 
and h/D as the abscissa. Thesecurves are shown.in Fig. 6. The hoop 
stress failure pressure was also plotted on these curves using a nominal 
value of compressive strength of 70,000 psi as obtained from DTMB Report 


1413.°° 


These curves were then used to select appropriate values for 

L/D and h/D. It was decided to keep the value of L/D constant at about 
2.5 and vary the value of h/D from about 0.05 to 0.09. It was also 
decided to test a minimum of three models of each geometry in order to 
determine that test results could be reproduced. Table 5 lists the models 
received from Zenith and gives their measured average dimensions, Fig. 7 
is a sketch of the model showing the end closure plugs. The plugs, made 
from medium steel stock, were designed to give simple support of the model 
ends, the condition for the von Mises buckling equation. The effective 


length, L, is as shown. The inside diameter was chosen as 3.195" because 


of the availability at the Zenith plant of a mandrel of that diameter. 
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Model Fabrication 

The models were built up using layers of glass fiber reinforced 
epoxy resin tape. The glass fibers ran in one direction, along the 
length of the tape strip, and the tape was applied to the mandrel both 
by circumferential winding under tension (circumferential plies) and 
by laying on strips parallel to the axis of the mandrel (longitudinal 
plies). The process was started and ended with two circumferential plies 
and was continued using alternate layers in a ratio of two circumferential 
plies to one longitudinal ply. This 2:1 ratio of circumferential to long- 
itudinal plies results in the allowable stresses in the longitudinal 
and ciremferential directions being in close agreement with the induced 
stresses in these directions. The resulting structure is of course aniso- 
tropic. The thickness was varied by varying the number of layers of tape. 
Five tubes were fabricated, each about 60" in length and each of a different 
thickness (see Table 5). Upon completion of the curing process, most of 
the models were cut from the tube to proper length while still on the 
mandrel at the plant. Some of the models were cut from the remaining 
portions of tubes after delivery, using a band-saw, and grinding or 
machining the ends until smooth and normal to the cylinder axis. By 
weighing and measuring some of the models, the density of the material 
was determined to be about 0.07 pounds/cubic inch, Values of resin 
content for the model material were provided by Zenith, and are given 


in Table il. 
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FIGURE 7 
































TABLE 5 
a 
3.364 
. Reg 
3.605 3.400 0.205 7.04 3.54 O06027 2 ay 
| 7.00 “ 6.50 . | 1 9f 
3.675 3.435 0.240 7.17 8.67 0.06987 2.52 
5.50 5.00 1.46 | 
3.757 3.476 O28/ 7.32 8.82 0.08084. 2S6. | 
4.62 4/2 L619 | 
3,813 3.504 0.309 9.44 8.74 0.08818 
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ALL DIMENSIONS ARE /{N INCHES 








Instrumentation of Models 

Strain data were taken to provide a means for calculating an 
effective modulus of elasticity for the models. As shown in Appendix C, 
values of the strain-pressure ratios, taken in the directions of principal 
stresses, can be used in conjunction with Hooke's Law and thick shell 
theory to solve for Eg and E,. A value for Poisson's ratio must be 
known or assumed. 

In order to determine average values for €,g and ¢€,, strain gages 
were applied to a number of the models, oriented both circunferentially 
and longitudinally. The number of gages used was limited both by 
financial considerations and by the number of leads that could be fed 
through the pressure fittings in the tank. In general, four gages were 
applied to each instrumented model. These gages were located at the mid- 
length of the model to avoid end effects, and were oriented two in a 
circunferential direction and two in a longitudinal direction. The two 
circumferential gages were placed "back to back", inside and outside, as 
were the two longitudinal gages, 

The choice of gage types was governed by: 

(1) the need to have non-pressure-sensitive gages; 

(2) the necessity for measuring strains up to 15,000 microinches/inch; 


(3) the desire to use a gage length long enough to reduce the effect 
of local variations in strain; and 


(4) the cost of gages. 

On the basis of these considerations, two types of gages were selected. 
For use on the inside of the model, where the gage would not be subject to 
the pressure of the fluid, the SR-4, A-5-1 wire page was used. This gage 


has a 1/2" gage leneth, is relatively inexpensive, and is suitable for the 
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measurement of strains up to 14 to 2%. For use on the outside of the 
model, where the gage would be subject to external pressure of the fluid, 
the SRe4, paper-backed constantan foil gage FAP-50-12 was selected. This 
gage, with a 1/2" gage length, is not as pressure-sensitive as the A-5-] 
gage. It also is suitable for measurements of strain up to 14 to 2%, 

Compensating gages of the same type and lot number were used in con- 
junction with the active gages. 

The gages were applied using Duco cement, and the gage manufacturer's 
instructions were followed. After the cement had dried, all gages were 
connected in the strain indicator circuit and pressed with a pencil eraser 
tip while watching the strain indicator to check for excessive deflections 


of the needle. 





Test Procedure 

Strain gages were applied to the inside of the rodel by using a 
"harness"=\{_a sheet of graph paper to which the gages were attached. This 
sheet was trimmed to exact dimensions, rolled into a cylinder, and positioned 
inside the model with glue applied to the gages. Leads of #28 stranded 
wire were attached to the strain gages. 

The ends of the model were prepared by thorough cleaning with acetone. 
PR-132]1 Class A sealant * was then applied to the ends, and the end plugs 
inserted, leading the wires from the inside of the mode] out through the 
1/4" tube ** screwed into the lower end plug. After the sealant set up, 
the model was ready for installation in the tank. Fig. 8 shows a model 
ready for installation and Fig. 5 shows a model in the tank and partially 
lowered into testing position. Fig. 3 shows a model*inside the chamber, 

The model was positioned in the open tank by first leading the wires 
from the 1/4" end plug tube through the "A" fitting in the center of the 
tank bottom, then connecting the external strain gage leads to the leads 
inside the tank from the "B" fitting. The model was then lowered, allowing 
the 1/4" tube to slide through the "A" fitting until the lower end plug 
rested on a steel "spool" support about 2" high located in the bottom of 
the chamber. The lock nut of the "A" fitting was then tightened, taking 


eare not to over-tighten. " 


* Manufactured by the Products Research Co., 410 Jersey Ave., Gloucester 
City, New Jersey : 


¥* This tube is a semi-permanent attachment to the end plug, its threaded 
end having been coated with sealant prior to screwing into the tapped 
end plug. 


§ 
f 


#*## Tt was observed that excessive tightening of the lock nut permanently 
deformed the ferrule and tube, making it impossible to slide the tube 
back out of the fitting after the test. 
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With the model in place, the strain gage leads coming from the 1/4" 
tube and from the "B" fitting underneath the pressure vessel were connected 
to the strain indicator or switching unit. The chamber was filled with trans- 
former oil * to a level slightly above the sealing surface. The circular 
plate was inserted, causing the excess oil to be squeezed into the plug 
recess, This precaution was taken to insure that the oil was relatively 
air-free, The screw plug was then positioned and tightened. Tapping the 
screw plug wrench handle with a hammer insured firm seating of the circular 
plate. 

Pressure was raised by operating the hand pump. The valve "G" in the 
line to the 5000 psi gage was shut at 4500 psi. This valve, although rated 
at 3000 psi, performed satisfactorily. 

In some cases, only one strain gage was measured per run, using the 
strain indicator directly. In most of these tests, at least two runs were | 
made per gage, and the readings for the two runs were found to vary very 
little. For other tests, the switching unit was employed and a number of 
gages were read per run. ) 

Upon completion of a test, the screw plug was removed, and the "A" 
fitting lock nut was loosened to permit removal of the model. As the model 
was removed, the 1/4" tube slipped out of the "A" fitting, allowing the oil 
in the tank to drain into a receptacle below. | 

The two pressure gages used during the tests were calibrated at the 
Material Laboratory before beginning the tests and at DTMB toward the end 
of testing. Both calibrations were made using dead-weight gage 


ee eee ee wee ee eee eee eee eee eee eee ees ees ee ese ee 


* WEMCO Class C transformer oi] was used as the pressurizing fluid, making 
it unnecessary to waterproof strain gages exposed to the fluid. 
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C= 





testers. No differences were noted in the gage calibrations. The 
5000 psi gage was adjusted after the second calibration to give gage 


readings closer to the actual pressure. The gage calibration curves 


are shown in Figures 9 and 10. 
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EXPERIMENTAL DATA 


A total of twenty-three models was tested, twenty-one of which were 
tested to failure. Seven of the models were instrumented to measure 
strains in order to determine the moduli of elasticity. Three models 
were designed to fail by yielding of the material. This was done in order 
to calculate the maximum compressive stress which the material would 
sustain and to determine the appearance of models after yield failure. 

Since the calculated circumferential stress for the "D" and "E" models 
was close to the nominal strength of the material, this information was 
used to ascertain whether these models did, in fact, fail by instability 
@lone, or whether some yielding might have taken place. 

A summary of the models tested, with the respective collapse pressures, 
is presented in Table 6, The recorded strain data are in Appendix A, and 
the pressure~strain plots are in Appendix B. Table 7 summarizes the 
measured strain-pressure ratios as obtained from the plots. The first 
models instrumented were tested with the compensating gage located inside the 
model for inside active gages, and with the compensating gage located in the 
pressurized fluid for external active gages. As shown in Appendix C, it is 
theoretically possible to cancel the effect of radial strain on the outside 
active gages by having the compensating gage mounted on a block which is 
subject to the hydrostatic pressure, Furthermore, any effect of pressure on 
the active gage in the fluid would be comrensated for , as well as more 
accurate temperature compensation provided. The later models tested, however, 
had all gages--inside and outside--balanced against a single dummy gage in- 
side the model. This latter procedure was followed because it had been 
determined that the error caused by neglecting the Poisson effect of the 


radial strain on the outside gages was negligible, and therefore the 


additional complication to the testing procedure was unjustified. 


IG 


In order to be able to compare and average the strain data taken under 
the two different methods of compensation, the strain-pressure ratios of 
active gages read with the dummy under pressure were corrected to values 
equivalent to those taken with the dummy gage inside the model. This was 
accomplished by simply adding the value of strainepressure ratio of a gage 
mounted on a dummy block under pressure, balanced against a dummy gage inside 
a model, All the strain-pressure ratios used in computations, therefore, 
have values which would result from active gages being read against a 


dummy located on the inside of the model, not under pressure, 








TABLE 6 
LAME STRESS IN MODELS AT COLLAPSE PRESSURE 
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MEASURED STRAIN- PRESSURE RATIOS 
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ANALYSIS OF EXPERIMENTAL DATA 


Strain Data Analysis 


The moduli of elasticity were calculated from the strain-pressure 
ratios using the formas presented in Table 8. The theoretical basis 
for these formulas is shown in Appendix C. A summary for the calcula- 
tions for Ey and E, is given in Tables 9 and 10. No values were cal- 
culated for the B-model because premature failure of the instrumented 
model] in this series prevented the taking of sufficient strain data, 
The values of Ey and E, used for the B-models are averages of the values 
for the other models. 

Since available instability formulas apply to isotropic materials, 
it was necessary to develop a method for converting the principal mod- 
uli Ey and E, into an effective modulus Eee This was accomplished by 
first deriving an expression for the shear modulus of an orthotropic 
material, The derivation followed a method used by Timoshenko? ° (pp. 
54-57) in deriving the expression for the isotropic case and is presen- 


ted in Appendix D. The shear modulus thus found is 


G= EXE g 
(Ex + Ey)(1 +7) 


Note that for E, = Eg this expression reduces to the isotropic case 


= E 
2(1 +7) 


In the derivation, Poisson's ratio was assumed to be constant. As noted 
later, this assumption does not lead to an appreciable error. Then, to 
obtain an effective modulus, the orthotropic expression for G was equated 


to that for the isotropic case as follows; 


a 











E 
Let Ex = kE,, then G = ~———__*_—-_ ,_ Equating this’ to 





1+k 
l1+y 
2 ( ) 
a Sg gives 
2(1 +7) 
k 
= 2B 4 eee 
Be 7 Sea 


The calculated effective values for the modulus are summarized in Table ll. 
In addition, the authors undertook the derivation of an expression 
for the buckling pressure of an orthotropic cylinder. This derivation is 
shown in Appendix E, and the buckling pressure is given in the form of a 
determinant. The calculation of buckling pressure using this determinant 
ig very cumbersom. Since the expression was derived shortly before the 
completion of the thesis, no evaluation of the determinant was made with 
respect to the experimental results contained herein. It is believed, 
however, that it may have merit in computing buckling pressures for ortho- 
tropic cylinders. Also, it might be used to verify the assumption of an 
effective modulus for use with the isotropic equations. It is noted that 
the determinant derived for the orthotropic case reduces to the expression 
for the infinite tube, with E, the only modulus, when L is made to approach 
infinity. This suggests that, in reality, an effective modulus for use in 
the buckling equation should depend on the length of the cylinder. It is 


also noted that the determinant reduces to that given by Timoshenko?4 
E 


(p. 449) when E, is made equal to E, and G is put equal to ~————_ 
2(1 +z) 


v 


except that the term - Pe ad is added to the last term in row 3, 
column C to take care of the end thrust. 

Details of the assumptions employed in the calculation of an effec- 
tive modulus follow. 


(1) As previously stated, Poisson's ratio was assumed to be constant. 


ae 





(2) 


e) 


On the basis of data provided by the David Taylor Model Basin 
and the Zenith Plastics Company, this value was assumed to be 
0.15. The effect of the variation in Poisson's ratio on the 
value of the calculated modulus is small. For example, using 
7 = 0.10 instead of 0.15 in the calculation for Eg raises the 
calculated value about 3% for both the A models and the E models. 
The stress field in the model was assumed to be in accordance 
with the Lamé thick shell theory, and the strain-pressure 
curves were assumed to be linear up to the point of collapse. 
As a consequence of this assumption the bending stress due to 
the simply supported tube of finite length is neglected. As 
will be shown later, the bending stress prior to buckling is 
a@ very small quantity. 

The effect of radial strain was neglected in calculating Eg, 
and Ex. This effect is small as is shown below: 


The radial strain on the outer surface is given by: (see Appendix C) 


ei, = Ee ee ee) 


Thus, 
{ K,/2 Ke 
IM. = Se eee ee 
2 ar ex re) 


It is seen that the greatest error would occur for the largest 
value of my, which value would correspond to the thickest 
cylinder, i.e., the cylinder with the least value of K). For 
the E-model, K, = 6.5; Ky = Ky - 1 25.5; By = 5.5 x 10°; 
By eae 10°; V= 0.15; then assuming a nominal value for 
E-=3x 10°, ec 1/3 - (0.15)( 2.28 + 235) = 0.046 


Now the exact formula for computing E using outside strain 


K2 (t+2¥)(i-2v) 
me, (1-Y¥)+Umy, +UMe2, 


gage data is: = 
Ea 





It was found that the difference between the value of Eg 
calculated by using this equation with m. = 0.046 and 
that calculated by neglecting m. entirely was less than 1%. 

It was further noted that the longitudinal compressive strain was 
greater on the outside than on the inside of the models. This difference 
is attributed to the fact that some bending of the shell undoubtedly 
occurs even for the lower pressures, That this difference is not inherent 
in the strain relations can be shown by the following development: 


From Appendix C, 


ae 7 
Ske 5 z( ee - ee (radial stress = 7 ) 


E a 1M 2S) 


hae Ey (radial stress = 0 ) 


Since Kp = Kj-1l, 


Z _ PK: _ were ve _ ver. wk ~upke 
Ex, — Ex, = + 2Fe a ie 


2 Ex ZEx Eg 
€ — € = oe Papeete [e: 
Re) Xs Eg Er 
from which, on 
/ 
My, - My, = “VEg =) 
Since E, =5x 10° E, =3x 10° and Y= 0,15, 


m,, - Mx, = (0.15)(1/5 + 1/3) = -0.02 
This difference is enall in comparison with the differences actually 
measured and is of opposite sign. It is concluded, therefore, that the 
measured difference results from bending of the generators. In order to 
reduce the errors arising from this effect, an average of m,, and mx, 


was used in all calculations which called for either value. 
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Table 7 shows that even for identical models there were considerable 
differences in measured values of strain-pressure ratios. These differences 
may be attributed to two causes: 

(1) The properties of the individual models may vary slightly, even 
though cut from the same tube. 

(2) Even on the same model, the strain-pressure ratios will vary from 
point to point around the circumference, as shown in reference (38). 

The most accurate data for use in computing moduli would be that taken 
from a Pinter of strain gages spaced around the circumference of each 
model, using the average values of strain-pressure ratios, Due to limi- 
tations discussed previously, this procedure could not be followed. 

Instead, the averaging process was accomplished by using inside and outside 
Strain data, and averaging the moduli thus calculated. Where circumferential 
variations in strain-pressure ratios are caused by the tendency to lobe 
formation, the use of inside and outside data from “back-to-—back" gages 

_ tends to cance] the variation. In the case of the A, D, and E geometries, 
two instrumented models of each were tested in order to provide more 

data for averaging. 

Since the tubes were fabricated using the same techniques and material, 
one would expect the moduli to be about the same for all the models. Such 
is the case, except that the A models were found to have a slightly higher 
modulus.(See Tables 10 and 11). As previously noted, the effective modulus 
for the B models was taken as the average of that of the other models, due 
to insufficient strain data. The small differences in resin content for 
the models (Table 11) are not sufficient to account for the differences 
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TABLE 8 
SUMMARY OF THICK SHELL FORMULAS FOR CALCULATING 


MODULI OF ELASTICITY 


Inside gages, dummy block not subject to pressure: 


aa : K,(1-v4 


a= fa 
? a Z2(my,+2/Mg,) 


Mg + 2”) »x, 


Outside gages, dummy block not subject to pressure: 


K2(1+%)(1-2y) Kk, (it vys-w 


EE = i a a 
a Ng, (1-2) ru(nx,.+M,) 2[mx2(1-v)+ v(m. + My)] 


Outside gages, dummy block subject to hydrostatic pressure**; 


7 (K2-1)1-v*) c, ae) 


Eg a / / / / 
Mb. +P Myer My2 + YUM¢g, 


## Dummy block assumed to be of same material as model; circumferential 
active gages balanced against circumferential dummy and longitudinal 


active gages balanced against longitudinal dummy. 
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TABLE 8 (Continued) 


SUMMARY OF THICK SHELL FORMULAS FOR CALCULATING 


Kae 


11g, 


/ 
Nx 2 


" 


MODULI OF ELASTICITY 


Circumferential modulus of elasticity 
Longitudinal modulus of elasticity 


Mean radius, '/2(R, +R.) 


Wall thickness, (¢.-R,) 


2R3 
Ri- R? 
RitRe- 
(mn eae 
Avg. strain-pressure ratio, inside circumferential gages, 
dummy not under pressure. 


Avg. strain-pressure ratio, inside longitudinal gages, 
dummy not under pressure. 


Avg. strain-pressure ratio, outside circumferential gages, 
dummy not under pressure. 


Avg. strain-pressure ratio, outside longitudinal gages, 
dummy not under pressure. 


Avg. strain-pressure ratio, radial, on outer surface. 


Avg. strain-pressure ratio, outside circumferential gages, 
dummy under hydrostatic pressure. 


Avg. strain-pressure ratio, outside longitudinal gages, 
dummy under hydrostatic pressure, 


Poisson's ratio. 
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TABLE 9 


DIMENSIONS OF TUBES FROM WHICH MODELS WERE TAKEN 


AND CALCULATION FOR K, 
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TABLE Il 


CALCULATION OF EEFECTIVE 
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Stress Analysis 


In order to properly evaluate experimental results, it was necessary 
to verify the mode of failure for each model. Although all models were 
designed to fail by instability, the predicted instability failure pressure 
for the D and E models was fairly close to the corresponding compressive 
yield failure pressure predicted using the thick-shel] theory and a nom- 
inal value for the yield strength of the material. Also, the thick-shell 
formulas used to calculate the stress do not take into account the bend- 
ing effects. Because of these considerations, the following steps were 
taken in order to make an evaluation for the mode of failure. 

(1) In order to evaluate the effects of bending, a formula was de- 
rived to calculate the longitudinal bending stress for a simply 
supported cylinder. The method of derivation was suggested by 
Mr. John Pulos of DTMB and is presented in Appendix F. The ex- 


pression for the bending stress at midbay is: 


Ea Bs [en* 7) ] (2) = na) es , S/N O42 
. (/- #4 L ie 


x=0 
Substituting values for the F model, the theoretical longitud- 
inal bending stress at midbay was found to be only 75 psi. 

(2) The Lamé (thick-shell) stress corresponding to the failure pres- 
sure was calculated for each of the models and compared both 
with a nominal compressive strength”of 70,000 psi as obtained 
from DTMB Report 1413, and with the lowest calculated value 
obtained from a yield model (63,233 psi). For the E model, the 
highest calculated stress was 59,285 psi (see Table 6). This 
value is 10,715 psi below the value for nominal strength but 


only 3,948 psi below that obtained from the yield model. It 


* Yield strength in circumferential direction, gg. 
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CosH 204, COS* 87, +. 51N 407, SIN*O ‘| : 





should be noted, however, that application of the Hencky-von 
Mises failure criterion would cause the differences to be greater. 
(3) Finally, the physical appearance of the collapsed models was 
compared to that of the thinnest model, for which a pure in- 
stability failure was certain, and with that of a short model 
which definitely failed by yielding of the material. Figures 11 
and 12 show collapsed C, D, and E models together with a short 
model which failed by yielding. It is noted that the longer 
models pictured failed along a longitudinal line parallel to 
the model axis by buckling inward, whereas the short model tended 
toward failure throughout most of its wall area, 
On the basis of these considerations it was determined that the A, 
B, and C models failed purely in the instability mode. Since the com- 
pressive stress for the E models was within about 4000 psi of the value 
obtained from a yield model, there still may be reason to question whether 
or not some yielding took place. Also, there was scatter in the points 
for the D and E models (See Fig. 13) which may be a result of local yielding 
caused by stress concentrations due to the heterogeneous nature of the 


material e 
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PRESENTATION AND DISCUSSION OF RESULTS 
IO LUN OF RESULTS 


Results of the experiments are shown graphically in Figure 13. The 
critical pressures determined from tests were converted to the non-dimensional] 
form p/E, . These values vere plotted against the thickness—diameter 
ratio, Calculations for p/E, are shown in Table a2; 

Also plotted on Figure 13 are curves representing various theoretical] 
buckling formulas, using an average value of L/D for the models tested, 
(Average L/D = 26D hs 

In spite of the anisotropic nature of the model material and the 
inaccuracies inherent in the computation for an effective modulus, very 
good agreement with theory was obtained exceptfor the thicker models. The 
fact that the points for the D and E models fall below the curves may indicate 
that the thin shell buckling theory is not adequate for ratios of h/pD in 
this range. As previously noted, however, some yielding may have occurred 
in these models, 

It appears that the von Mises equation is the best one to use in the 
range of h/D> 0.05 Since, in general, it yields conservative values for 
the collapse pressure, Design curves for the buckling of simply supported 
cylinders were therefore draw up for the von Mises equation, (See Fig, 14). 
Data for these curves were calculated using the IBM 7090 computer at DTMB. 
These data are reproduced in Appendix G, 

Since curves were not available for determining the number of circum- 
ferential lobes into which a relatively thick cylinder would fail, curves 
were drawn. (See Fig, 15). These curves were calculated using data in 
Appendix G by putting the buckling equation in the form 

P/E = Fy(h/D)3 + Fo(h/p) 
for each n and L/D and then solving for the desired intersections, These 


Curves are for y= Clore 


eee Oe 









For values of n with Y = 0.3, in the range h/D€ 0.02, see reference 
(33). 

As far as is known by the authors, these curves for the von Mises 

equation and for the number of lobes are the only ones available in the 


range of h/D greater than 0.05. 
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CONCLUSIONS 


Based on the results of these experiments the following conclusions 


are drawn: 


(1) 


(2) 


(3) 


(4) 


The von Mises equation, equation (6), Ref. 22 is the best in- 
stability equation for cylinders with ratios of h/D greater than 
0.05. It gives conservative values for the buckling pressure up 
to an h/D of about 0.07. Use of presently available simplified 
forms of buckling equations cannot be expected to give conserva- 
tive values for the buckling pressure in this range of h/D. It 
is noted that Windenburg and Trilling’? had previously reached 
the conclusion that the von Mises formula was probably the best 
one; however, their conclusion was based on an analysis of exper- 
imental and theoretical results for cylinders with ratios of h/D 
less than 0.007 in contrast with the much thicker range of 0.05 
to 0.088 investigated in this thesis. 

For practical cases, the use of stiffening rings, heavy webs, and 
bulkheads would result in smaller ratios of h/D and L/D for a 
given weight of structural material. It is concluded, then, that for 
most practical cases, the thin shell theory is sufficiently 
accurate for predicting instability failures. 

For cylinders built up of fiber glass reinforced plastic layers, 
the methods employed in this thesis are sufficiently accurate 

to predict instability failure pressures. 

For an orthotropic cylinder, the effective modulus, as used in 
the buckling equation, is probably a function of the length 


between supports. 
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RECOMMENDATIONS 


Further Investigations, 
The experimental apparatus designed and assembled for this thesis 


provides a means for the undertaking of a variety of experimental in- 


vestigations in the field of hydrostatically loaded shells. Some investi- 


gations which might be undertaken are; 


(1) 


(2) 


Investigate the stresses in short, relatively thick, fiber- 
glass reinforced cylindrical shells. Measure the strain 
distribution along a generator and determine the effects of 
bending. Verify the formula for bending stress which is 
derived in Appendix F, 

Using fiber-glass reinforced plastic models, verify the ortho~ 
tropic buckling determinant as derived in Appendix E. In 
conjunction with this, evaluate the determinant to arrive at 
a general expression for buckling pressure similar to the 

von Mises equation for the isotropic case. Using the result, 
develop a method of computing an "effective modulus" for use in 


isotropic buckling formulas. 


(3) Perform experiments on models of unusual form, such as: 


a. Stiffened spheres 
b. Ellipsoids 


ec. Sphere-cylinder intersections 
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RECOMMENDATIONS (Continued) 


Improvement of Apparatus, 
Replace valve (G), (Fig. 4) with one rated at 10,000 psi or higher, 


such as American Instrument Company valve, Cat. No. 44-1505. 

Install a fitting in the line to each gage to damp out the sudden 
drop in pressure when a model breaks, thus protecting the gage pointer. 
There are commercial fittings available for this purpose. 

If a greater number of strain gage leads are required to be run for 
the gages mounted on the outside of the model, it is recommended that a 
CONAX Corporation thermocouple fitting, Cat. No. TG-20-B8, be installed 
in place of the fitting "A" now installed. The TG-20-B8 will take eight 


#20 wires. 
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APPENDIX A 


The following tables consist of the recorded strains corresponding 
to the gage pressures and corrected pressures. The model number, gage type, 
gage factor, date of test, gage location , cycle number, run number, and 
pertinent remarks are also indicated. 


Three types of Baldwin-Lima-Hamilton gages were used for the strain 


measurements: 
Gage Type Nominal Gage Length 
A=-5-1 1/2" 
FAP=50-12 1/2" 
FA~100-12 i 


The type A-5-1 gage was used on the inside of the model only, since 
the wire gage is particularly sensitive to hydrostatic pressure. The 
FAP-50-12, a foil gage, was used both inside and outside the models. The 
FA-100-12 was used on the outside of model A-4 only. All three types 
appeared to give reliable data. No particular difficulty was experienced 
in handling and applying the gages. 

In the row entitled Gage Location, "L" stands for a longitudinally 
mounted gage, "C" for a circumferentially mounted gage, "1" for an inside 
gage, and "2" for an outside gage. 

In some cases the gages were read one at a time requiring a differ- 
ent pressure cycle for each gage. A record was kept of the total number 
of cycles for each model and the number of runs for a particular gage. 
These data were taken in order to ascertain the effects of creep, if any. 
It appeared from the results, however, that creep was not an important 


factor for these short term tests. 
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APPENDIX B 
PLOTS OF PRESSURE VERSUS STRAIN 


This section contains plots of the measured strains. The resulting 
Slopes, designated by the letter m with appropriate subscripts are indi- 
cated on the plots. 

It is noted that all plots are very nearly linear until the col- 
lapse pressure is approached. The reversed slope for gage L-2, model 
A-4 indicates the development of a tensile stress as a result of bend- 


ing near the collapse pressure. 
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APPENDIX _C 
THEORY FOR CALCULATION OF Ey AND Ex FROM STRAIN DATA 


The circumferential stress in a thick-walled cylinder under external 


hydrostatic pressure, p, is given by: 


2 Ricken 
Fe eos a a L rer (42) | 


To THE POINT 


where r is the distance from the center, at which the stress is measured, 
If r is given values of r = Ry and r = Ro, the following expressions 
evolve for the circumferential stress on the inner and outer walls of 


the cylinder; 


2Ri 
Co ee pEEe 


| 


es 2 Rs 
io. = PCRar) = Paar ~!) 


(gd, and Oz, are positive if compressive. ) 
The longitudinal stress is constant over the cross section, and 


is given by; 


Re 
= ? Ri_-R* (Ref. 42) 


aes 


——— ») we have 
A -R, 


ii 


Letting kK, 


Og, = OK, 
Tb p(K,-1) = pXa (1) 


\\ 


If we assume a material which has three principal values of Young's 


Modulus and a single value of 2, then from Hooke's Law: 


= Ge 





ow eee ee 











ena Ox ee a 2.20, (2) 
Ex ae a, 
Ox ZOK  _ wg 








ae are Eg 


Gombining Hooke's Law and thick shell theory, we have strains in the . 
inside and outside surfaces of a thick-walled cylinder Subject to external 


hydrostatic pressure: 











» Kame eR) ‘ 
ae rE [et Ee 
€x2 = pa a aa (OUTSIDE ) (3) 


a f Be) 
Er2 - Ge ZF ES 





(Note that c= p on the outside surface) 














ion a) 
Eg, = 1 Eg@ oe 2£Ex 
Ki t 
ex, me P wee ae (iv SIDE) (4) 
K, iy 
oS See Ea? 


(Note that o,= 0 on the inside surface) 


The stresses in the surfaces of a block of material subject to 
external pressure, p, are given by: 
og = oO. = Oo, = Pp 
From Hooke's Law, the strains on the outside surface of the block 
are given by: 
| ya ees ; 
Egp = ia Erp Exp 


= Oe (S) 





Y Lv 
se Erp Exp rs 


The > subscript on a term indicates that it refers to the block rather 


than to the cylinder. 


eee 





Consider the case of strain measurements made with strain gages 
on the inside surface of the cylinder, using a dummy gage not subject to 


pressure. The strains as read are given by: 


eee ae 
€¢, = phil 3E;) 
3 / Zo 
Cx, = pK GE - £5) (G) 
ae A(t Bae 
er = PUK. (zee ea) 
Solving the first two of these equations: 


Mae (1 —z7*) 
Mg +2 Mx, 


Eg = 
(7) 


re) 
a (Ny, +VIM¢,) 


fF. = 


Where iy Se 2 


/ 


a 

Therefore a model instrumented thus can be subjected to hydrostatic 
pressure and the measured strain~pressure ratios, mand mg, , with 
an assumed or known value of z can be substituted to solve for E, and Ey. 
The radial strain does not enter the solution, and the radial stress is 
zero. 

Consider the case of strain measurements made with strain gages 
on the outside surface of the cylinder, using a dummy gage not subject 
to pressure, The strains as read are given by: 

és = Pp SS (8) 


ZE* aoe ae 
= (2 eee 
a ae 7 











ea 


Solving this system of equations simultaneously gives: 


-C3- 








fp = eA 


Co [ae arene 
re 
pe 


ae Se ee $2+Mr,)| 


Consider the case of strain measurements made with strain gages 
on the outside surface of the cylinder, using dummy gages mounted on a 
block subject to the hydrostatic pressure, the longitudinal active gage 
being balanced against a longitudinally oriented dummy, and the circum- 
ferential active gage being balanced against a circumferentially oriented 
dummy. Assume the dummy block to be of the same material as the cylinder, 


The indicated strain readings are; 


A 
Co. = Cg2 — €¢p 
Ex, = Exe —€E¢gp 
/ : 
See ira ee 


where the primes indicate that the strains are the indicated values rather 


than the actual strains. From (3) and (5): 














= ES 2 i eee ay, | aa (aye 
By2 F E¢ ee eS Eur =, 
ee of eee ee ee oem en 

e725. Ee Ee / vss Er Ey, 


f 


a > —— = ee P eae Re pe 


(Az-t) _— o( Kx =) a 
4. = a ie. ) 
_- s -)_ o(hr-1) . » (10) 
x2 = a 
ax zs f 
/ 


aA = 





The terms involving radial strain drop out. A comparison of (10) 
with (3) shows that they are identical, except that in (10) the radial 
terms are missing, and Kp - 1 and K,/2 - 1 appear in place of Ko, and 
K,/2, respectively. Thus by having the compensating gages arranged as 
in this case, the values of Eg and ES may be computed from the indicated 
strains €4, and Exes and the radial strains do not enter the solution. 


Solving (10): 


(Ke-1)(i-2v*) 
Mg, + UM xy 


_ (K./2 -I)C1-v) 


Mx, + UMB, 


£¢ 
(11) 


Ex 


WHERE : Mg, = €p2/ 


Mya = €x2/v0 


-c5§ - 











APPENDIX D 


DERIVATION OF SHEAR MODULUS FOR AN ORTHOTROPIC MATERIAL 
Following a method similar to that used by Timoshenko (Ref. 40, 
pp. 54-57) for developing the shear modulus for the isotripic case, a 
formula is derived for the shear modulus of an orthotropic material 
by assuming two principal moduli of elasticity and a constant value 
for Poisson's ratio. 
Consider an element in pure shear, For this condition o% = - = Tike, 


as is readily apparent from Mohr's circle. The element can be represented 


as shown below: 
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(The deformed element is represented by dotted lines.) The deformed 


angle at A is3-Y; that at B is Fay » where Y is the shearing strain. 
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Solving for Y g/ves 
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NeglecTing The S@cond term ¢n the ASenonmsnaoJor and solving For 7xeo 
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Conway”! arrived at a similar expression using an energy analysis. 





APPENDIX E 


DERIVATION OF EQUATION FOR THE CRITICAL PRESSURE 
OF A SIMPLY SUPPORTED, ORTHOTROPIC, CIRCULAR 


CYLINDER SUBJECTED TO HYDROSTATIC PRESSURE 


Stress-strain relations. 
For isotropic shells, the relation between the stress resultants 


and couples and the corresponding strains is: 
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Nig = GY eo (O72) ~ xe 


where the subscript “o” refers to strains in the middle surface. 
If is is assumed that the cylinder is orthotropic with principal 
moduli Ex and Eg and with a constant value of Poisson's ratio, these 


equations may be written: 
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changes of curvature in terms of the displacements 
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The stress resultants and couples can then be written in terms of 


the displacements; 
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The Equilibrium Fquations. 
The equilibrium equations for a circular cylinder under radial 
load only are: 
2 2Nx 2 Ne 
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[see Timoshenko, Ref. 30, p.447; note that the term ga (arse x 


has been neglected, as was done by von Mises“*, | 
. These equations may be modified to include the effect of axial 
load as follows. This procedure was used by von Mises. 
The total axial force on the ends of the cylinder is equal to 77 a*p. 
This force is transmitted to the ends of the cylinder of area 27Qh. 


Hence the pressure in the cylinder walls is 
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Now consider an element of the shell in the longitudinal plane. 
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After deformation, the curvature of the generatrix (neglecting 


higher order terms) is 
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The force per unit length, pth, thus has a component in the z-direction 


of 
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Since the radius, a, had been divided out of the original equilibrium 


equation, this term is written as 
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In order to include the effect of the axial load, this term is added 


to the term in parenthesis in equation (3) to get: 
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Boundary conditions, 


The transverse and longitudinal cross sections of the cylinder assume 


the following shapes after buckling: 
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At x = 0, the displacements may be written: 
“=A snne 


na = B cosne wT 7s 90° ovt of Phase with « and «v7. 
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At x = 0, u = 0 and v and w are maximum; hence, the displacements for 


any value of x become; 
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Determinant for the Critical Pressure, 

By substituting the boundary conditions into the relations for the 
stress resultants and couples, and substituting the resulting expressions 
into the equilibrium equations, a system of three equations in A, B, and 
C is determined. In order that an instability condition exist, the constants 
A, B, and C, which determine the magnitudes of the displacements, must 
have values other than zero. This condition requires that the determinant 
of the coefficients of A, B, and C be zero. The determinant for this case 
ig given on the next page. Evaluation of the determinant for a given 


geometry and material yields the critical buckling pressure. 
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APPENDIX F 


DERIVATION OF AN EXPRESSION FOR THE LONGITUDINAL BENDING STRESS 
IN A 
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Salerno and Pulos43 have derived expressions for the stresses in 
a circular cylindrical shell supported by uniformly spaced circular 
rings of constant cross section. By using the displacement function 
obtained in their paper with the appropriate boundary conditions for 
simple support, an expression for the longitudinal bending stress at 
midbay is derived. 


The longitudinal bending stress is given by: 
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The displacement function,«u, as derived by Salerno and Pulos for 


deflection in the radial direction, is 
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where B and F are constants depending on the boundary conditions, A, # A3, and 
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Taking the x-axis along the axis of the cylinder, with the origin at 


the mid-length, the boundary conditions for simply~supported 
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Solving simultaneously, 
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Krenzke and Short ““ have presented a graphical method for calculating 


the stresses obtained by the analysis of Salerno and Pulos. 


In this 


reference, different parameters were employed as a simplification. 


These parameters, defined below, are used for the remainder of this 
derivation, 
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Using the Krenzke and Short parameters, 
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Taking the second derivative of this expression, simplifying, and sub- 


stituting into the expression for the longitudinal bending stress gives: 
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APPENDIX G 


SOLUTION OF VON MISES' EQUATION (6) 





This appendix contains solutions of the von Mises equation (6)* 
as computed by the IBM 7090 computer at DTMB. Values of p/E were com- 
puted for a range of L/D from 0.50 to 3.00 in increments of 0.25, and 
for a range of h/D from 0.01 to 0.10 in increments of 0.01. The solution 
for a particular h/D was started with a value of n = 2, The value of n 
was then incremented by one until a minimum value for p/E was obtained. 
Columns 3, 4, and 5 give the auxiliary quantities ale and 4, for 


use in the von Mises equation. 
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